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oo '■ Abstract 

It is an important question whether the final/initial state gluonic interactions which lead to naive- 
time-reversal-odd single-spin asymmetries and diffraction at leading twist can be associated in a 



Oh 



X 



definite way with the light-front wave function hadronic eigensolutions of QCD. We use light-front 
time-ordered perturbation theory to obtain augmented light-front wave functions which contain 



an imaginary phase which depends on the choice of advanced or retarded boundary condition 



for the gauge potential in light-cone gauge. We apply this formalism to the wave functions of 
the valence Fock states of nucleons and pions, and show how this illuminates the factorization 



properties of naive-time-reversal-odd transverse momentum dependent observables which arise from 
rescattering. In particular, one calculates the identical leading-twist Sivers function from the 
overlap of augmented light-front wavefunctions that one obtains from explicit calculations of the 
single-spin asymmetry in semi-inclusive deep inelastic lepton-polarized nucleon scattering where 
the required phases come from the final-state rescattering of the struck quark with the nucleon 
spectators. 

PACS numbers: 
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I. INTRODUCTION 



Wave functions are key objects of the quantum world, specifying the structure of com- 
posite states in terms of their fundamental constituents. The conceptual extension of the 
non-relativistic wavef unctions of Schrodinger theory to relativistic hadron physics are the 
frame- independent light-front wave functions (LFWFs) of hadrons ^^(xi, kij_, where 
Xi — -pip = pomp's are the light-front momentum fractions of the n constituents, ki± the 
transverse momentum components, and A, the parton helicities. The LFWFs are defined 
as constituent wave functions at fixed light-front time r = x + = t + zip and in the light- 
cone gauge A + = A + A z = where A 11 represents the gauge field l|-|4j. The LFWFs are 
obtained explicitly by computing the hadronic eigensolutions \^h) of the QCD light-front 
Hamiltonian H^p projected on the free Fock basis = (n\^fji). 

Light-front wave functions in QCD describe the quark and gluon composition of hadron 
at a fundamental level, leading to a description of a wide range of hadronic and nuclear 

n 

physics phenomena |2J. For example, the parton distribution functions measured in deep 
inelastic lepton-hadron scattering, including DGLAP evolution and their transverse mo- 
mentum extensions, are defined from the sum over squares of the light-front wave functions. 
Form factors are given by the sum of overlap matrix elements of the initial and final LFWFs 
with the electroweak currents. The gauge-invariant distribution amplitudes 0(2;, Q) which 
control hard exclusive reactions are the valence LFWFs integrated over transverse momenta 
k 2 ± < Q\ 

Recent theoretical developments have shown that final and/or initial-state interactions 
can generate a phase in scattering amplitudes which lead to novel single transverse spin 
asymmetries in high energy hadronic reactions at leading twist. A prime example of this 
rescattering physics in QCD is the Sivers single-spin asymmetry measured in semi-inclusive 
deep inelastic scattering and spin-dependent Drell-Yan lepton pair production fl-Q. Double 
initial-state interactions lead to an anomalous cos 20 azimuthal dependence of the production 
plane in unpolarized lepton pair hadroproduction, corresponding to the breakdown of the 
Lam- Tung relation in PQCD jsl, Similarly, diffractive deep inelastic lepton scattering 
£p — > i'p'X arises from the exchange of gluons in the final state which occurs after the hard 
lepton-quark interaction [loj]. Since nuclear shadowing involves diffractive deep inelastic 
processes, nuclear distributions are also dependent on rescattering processes. 
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The wavefunctions of stable hadrons that are obtained by solvin g th e Heisenberg problem 
H®p D \^ >= Ai 2 \^ > have a real phase. As discussed in Ref. 11], one can distinguish 
"static" structure functions, the probabilistic distributions computed from the square of 
the light-front wavefunctions of the target hadron from the "dynamic" structure functions 
measured in deep inelastic lepton-hadron scattering which include the effects of rescattering 
associated with the Wilson line. Thus it is an important question whether the final/initial 
state gluonic interactions responsible for the dynamics of rescattering can be associated 
in a definite way with the light-front wave function eigensolutions of QCD. The resulting 
augmented LFWFs provide an important tool for understanding the factorization properties 
of dynamical hadronic phenomena including single-spin asymmetries and diffraction. 

It has been shown that the light-cone gauge condition A + = does not fix the gauge of 
Abelian or non-Abelian gauge fields completely [7| : one has to choose a boundary condition 
for the transverse component of the gauge potential at spatial infinity: Aj_(x~ = ±00) 
The propagators of the gauge field which define the QCD Light-Front Hamiltonian in the 
Heisenberg problem are regulated using the principal value prescription. However, a different 
choice of boundary condition will lead to different properties of the light-front wave function 
amplitudes. In particular, if we choose a retarded (Ax(x~ = —00) = 0) or advanced 
(A±(x~ = 00) = 0) boundary condition, the resulting augmented light-front wave function 
will contain the necessary phase to generate the nonzero single spin asymmetry in hadronic 
reactions. We will demonstrate these properties, giving an explicit calculation in light- 
front time-ordered perturbation theory [21]. The result of our analysis provides the general 
structure of augmented LFWFs which is easy to apply to phenomenological applications. As 
an example, we will present results for the three-quark Fock state component of nucleon and 
the quark-antiquark component of pion at lowest non-trivial order. We can further simplify 
the result for the pion in terms of the distribution amplitudes. Given these light-front wave 
function amplitudes results, it is straightforward to calculate the pseudo-time-reversal-odd 
quark distributions of the nucleoli and pion. by applying the overlap formalism derived 
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121]. The light-cone gauge with retarded/advanced boundary condition has also been 



used to investigate the small- 2 physics 
factorization for nucleus-nucleus collisions 



13l4l6l| . in particular, to study the evolution and 



16|. 



The rest of this paper is organized as follows. In Sec. II, we present a general derivation 
of augmented LFWFs using light-front time-order perturbation theory within a lowest order 
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formalism. In Sec. Ill, we apply our method to the construction of augmented light- 
front wave function amplitudes for the three-quark Fock component of nucleon and quark- 
antiquark component of the pion. We summarize our paper in Sec. IV. 

II. GENERAL DERIVATIONS 

We start our derivation by constructing the general form for a Fock state expansion of 
any given hadron, 



where P and S are the momentum and spin of the hadron, d[i] = dxid 2 ki±/(y/2xi(2ir) 3 ) 
with the overall constraint on Xi and ki±_ implicit. For convenience, in following calculations 
we set the transverse momentum of hadron equal zero: P± = 0. Because the wave functions 
are boost invariant, all our results can be extended to more general case with P± ^ 0. In 
this Fock state, each parton is represented by the associated creation operator a\{ki) with 
hi = {kfik±i) = (xiP + , k±i), which contains certain longitudinal momentum kf = XiP + 
and transverse momentum whereas the minus component is determined by the on-shell 
condition k~ = {k\ + mf)/kf. Implicitly, the above light-front wave function amplitude ip n 
depends on the orbital angular momentum projection from the constituents with the form 
of (kf ± fcf) [171 ] . Since the following derivation does not depend on this structure, we will 
not include it explicitly. 

As discussed in the Introduction, the definition of a light-front wave function amplitude 
ip n can be extended to include rescattering effects so its phase is not necessary real. We 
can obtain the imaginary part (or the phase), by iterating the light-front wave function 
eigensolutions employing a particular boundary condition for the gauge field. In fact as 
we shall show, the phase of the augmented wavefunctions can be computed perturbatively 
by applying light-front time-ordered perturbation theory j2[ analogous to the Lippmann- 
Schwinger method. 

The first order correction to the LFWF can be obtained by iterating the Light-Front 
equation of motion: 




(1) 
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FIG. 1: Light-front time-order perturbation Feynman diagrams for the phase contribution from 
one-gluon exchange between two constituent quarks. 

where E ^~ represents the sum of all partons energy k~ , d[i]' represents the integral of 
{yiAii-)- The interaction kernel K can be calculated from the light-front time-order pertur- 
bation theory The wave functions tp n and ip' n may differ. From the above expression, 
we find that the phase of ipn m ay come from the wave function in the right hand side ip' n 



or the interaction kernel K. In the following, we assume that the wave 
for example, from model calculation such as constituent quark model 



unction ip' n is real, 



18] . We will focus on 



the contribution from the interaction kernel. We will calculate, in particular, the one-gluon 
exchange contribution to the interaction kernel. 

At the lowest order of the light-front time-order perturbation theory, we have one gluon 
exchange contribution to the interaction kernel. This can be expressed as a sum of all 
diagrams with gluon connection between all possible pair of constituents in the light-front 
wave function. For example, the contribution from the gluon exchange between the ith and 
jth quark can be written as, 



!h y/Xj v 

1 % 



P--q~-K-q- E k~ + ie <? 
1 6{q' + " 

p--q'--kj-er- E k- + *£ Q' + 



(3) 



where A represents the helicity for the associated quarks, q + = k~j~ — £~l and q' + = kf — 



and the color factors are implicit in the above equation. Similar expression shall hold for the 
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gluon constituent in the wave function, and so the final results. We illustrate the contribution 
in the above calculations for i — 1 and j = 2. The first term in the above bracket comes 
from Fig. 1(a), whereas the second term comes from Fig. 1(b). Moreover, at this particular 
order, quark number is conserved, such that n = n' . The gluon polarization tensor is defined 
as, 

dfa, = —g^u H [v • g] ' ^ 
where v is a light-like vector v ■ P = 1, q differs from q in the minus component to take into 

n 

account the instant propagator contribution |2J. When one solves the Heisenberg eigenvalue 
problem for light-front QCD as in discretized light-front quantization f3], the light-cone 
gauge singularity is regulated by the principal value prescription, which corresponds to the 
antisymmetric boundary condition for the gauge potential A±(x~ = +oo)+Aj_(x~ = — oo) = 
0. However, this prescription will not result into an imaginary part for the light-front wave 
function from the above interaction kernel. In the following calculation, we will choose the 
advanced boundary condition: Aj_(x~ = +oo) = whereas Aj_(x" = — oo) ^ in order 
to construct the augmented light-front wavefunction. With this boundary condition, the 
light-cone singularity will be regulated by [3], 

v^ + v^ v^q v v u q^ 

? 1 \Adv. = H — — , (oj 

[v ■ q\ v ■ q — ze v ■ q + le 

where the momentum flow of q is toward to the vertex v. Clearly, this term contains a phase. 

The imaginary part is simple, and proportional to a Delta function: m5(v ■ q). Since we are 

only interested in the imaginary part of the light-front wave function amplitudes, we simply 

apply this Delta function to the interaction kernel in Eq. (3). In particular, we find that the 

dominant contribution comes from the d + ± components of the tensor (3, [isj]. All other 

contributions cancel out between the above two terms or by themselves. Another important 

consequence is that the helicities are conserved in the interaction kernel: 5a;A' ■ After a 

little algebra, we obtain a rather simple result for the imaginary part of the light-front wave 

function amplitude generated from lowest order perturbation theory, 

i[M^k a± )] = -g[c.F.] (i - p"Z^fc- ) E^fo; V) . ( 6 ) 

where [C.F.] represents the color-factor for the Feynman diagram in Fig. 1 and ipn = 
if) n (x a ;£i± = ki± — q±,£j± = kj± + Q±,£p±\p^i,j = kp±). We emphasize that the wave 
function at the right hand side only contains real part. 



The above equation is the main result of this paper. It explicitly demonstrates that 
the light-front wave function amplitudes contain an imaginary part if we choose advanced 
boundary condition for the transverse component of the gauge potential. If we choose the 
retarded boundary condition, we obtain an opposite sign in the above equation. 



III. APPLICATIONS TO PION AND NUCLEON 



The three-quark Fock state components have been classified in Ref. 12]. For these light- 
front wave function amplitudes, we can apply the derivation in the last section, and obtain 
the imaginary part as, 

l[V qqq (xi,k±)] = ^C B J d 2 £ 1± d 2 £ 2± d 2 £ 3± 6 {2 \£ 1± + £ 2± + £ 3± )V qqq (x l ,£ l± ) 



p--5> 



5(2) ^-^ + (2^3) + (1^3) 



(k 1± -£ 1± ] 
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(7) 



where Cr = (N c + 1)/2N C and \I/ represents the general wave function amplitude constructed 



in Ref. 



121 ] . By applying these results, we are able to formulate the naive time-reversal-odd 



quark distributions (such as the quark Sivers function) in terms of the light-front wave 
function amplitudes, by taking into account the above imaginary phase using the above 
derivation [l2j. 

For the quark-antiquark Fock component of pion, the result can be further simplified as 



a Sn fd 2 q± l( , v A x(l-x)M 2 -(k ± -q ± ) 2 -m 2 g 

-C F J -g-^s, k x - g± ) [l - x{l _ x)M ,_ kl _ ml 



where we have chosen P± — and assumed that the quark and antiquark have the same 
mass m q and CV = (Nq — l)/2Nc- In particular, if we are interested in the large transverse 
momentum behavior of the light-front wave function amplitudes, we can expand the inter- 
action kernel in terms of £±/k±, and keep the leading order contribution. By doing that, we 
will obtain, 

ry 1 

I[^ x ,k ± )] = -^ 7 -C F( p(x) , (9) 

and 4>(x) is the leading-twist distribution amplitude for Pion, normalized by the leading 
Fock component light-front wave function <p(x) = J d 2 £±ip(x,£±). Similar expressions can 
be found for the quark-diquark model 0. 
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Although light-front wave functions depend on the boundary condition of the gauge 
potential in the light-cone gauge, physical observables cannot depend on this choice because 
of gauge invariance single-spin asymmetry in semi-inclusive deep 

inelastic polarized proton deep inelastic scattering £p^ — > i'qX and the associated quark 
Sivers function can be formulated simply as the overlap of augmented LFWFs using the 
advance boundary condition 12] . In particular, it is the phase difference between the LFWFs 
for the S and P-wave Fock components that contributes to the quark Sivers function in the 
quark-diquark model studied in Ref. The imaginary phases are calculated by using the 
general formalism Eq. (6) with similar expression as Eq. (8). 

The result for the Sivers single-spin asymmetry using augmented LFWFs is identical 
to that found in Ref. Q| using conventional LFWFs (with the principal value boundary 
condition), together with an explicit calculation of the final state phases which arise from the 
rescattering of the struck quark with the spectator diquark after the lepton-quark interaction. 
This identity is possible since the final-state phase due to rescattering is independent of the 
momentum transferred in the lepton-quark interaction. On the other hand, If we choose the 
retarded boundary condition, the augmented wave function will have opposite imaginary 
part. However, under this boundary condition, we have to take into account the final state 
interaction effects (the gauge link contributions from the quark distributions), but again, 
this leads to the same result compared to that using the advanced boundary condition. 

Similar conclusions hold for the small- a; parton distribution calculated in lOj. We leave 
this topic for a future publication. 



IV. SUMMARY AND DISCUSSIONS 



We have use light-front time-ordered perturbation theory to obtain augmented light- 
front wave functions which contain an imaginary phase which depends on the choice of 
advanced or retarded boundary condition for the gauge potential in light-cone gauge. We 
have applied these results to construct augmented wavefunctions for the three-quark or 
quark-diquark Fock state components of nucleon and the quark-antiquark component of the 
pion. We obtain the leading-twist quark Sivers function from these augmented light-front 
wavefunctions, by applying the overlap formalism [12!]. The result is identical to the explicit 
calculation [5] of the single spin-asymmetry in semi-inclusive deep inelastic lepton-polarized 



S 



nucleoli scattering where the required phases come from the final-state rescattering of the 
struck quark with the nucleon spectators. 
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